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\D • Abstract 

We report the results of a blind search for flavor-changing neutral current (FCNC), lepton- flavor 
violating, and lepton- number violating decays of D + , Df, and D° mesons (and their antiparticles) 
into 2-, 3-, and 4-body states including a lepton pair. Such decays may involve Flavor-Changing 
Neutral Currents, Leptoquarks, Horizontal Gauge Bosons, or Majorana Neutrinos. No evidence for 
any of these decays is found. Therefore, we present 90% confidence level branching-fraction upper 
limits, typically at the 10 -4 level. A total of 51 decay channels have been examined; 26 have not been 
O 1 previously reported and 18 are significant improvements over previous results. 



Introduction 

The E791 Collaboration has reported limits on rare and forbidden dilepton decays of charmed D 
mesons PJEHHj- Such measurements probe the SU(2)xU(l) Standard Model of electroweak interac- 
tions in search of new mediators and couplings [IJ [3] beyond the and Z° discovered at CERN in 
1983 pj- Here we summarize the results of two related analyses. First [2] we examined the tt££ and 
K£l decay modes of D + and Df and the £ + £~ decay modes of D°. Then we extended the methodology 
to 27 dilepton decay modes of the D° meson 3 containing either resonant V£ + £~ decays, where V 
is a p°, K* , or (ft, and non-resonant h\h?,£l decays, where hi is either a it or a K. The leptons were 
either muons or electrons. Charge-conjugate modes are implied throughout this paper. The modes are 
lepton flavor-violating (e.g., D + — > ir + fi + e~), or lepton number-violating (e.g., Df — > 7r~/i + ^ + ), or 
flavor-changing neutral current decays (e.g., D° — > K*°e + e~). WW box diagrams can mimic FCNC 
decays, but only at the 10~ 10 to 10 -9 level [Sl[7j. Long range effects through resonant modes (e.g., 
D° — > K*°p°, p° — > e + e~) can occur at the 10 -6 level Neither Lorentz nor gauge invariance [5] 
require lepton number conservation. A leptoquark ^H] m & n exchange diagram (e.g., D° — > p + e~) or 
a Horizontal Gauge Boson in a spectator diagram (e.g., D + — > tt + p + e~~) might mediate a change 
both in quark and lepton generation simultaneously. A Majorana Neutrino might lead to same sign 
dilepton decays [T^J. Numerous experiments have studied rare decays of charge -1/3 strange quarks. 
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Charge 2/3 charm quarks are interesting because they might couple differently |13| . 

The data come from measurements made with the E791 spectrometer |14j at Fermilab's Tagged 
Photon Lab. The spectrometer has been upgraded for a series of charm experiments including E516 
[T5| . E691 jTB], E769 ^7j> an d E791. E791 events were produced by a 500 GeV/c ir~ beam interacting 
in a fixed target consisting of five thin disks, one platinum and four diamond. The discs were well 
separated to allow charmed hadrons to decay in air spaces where other interactions would be minimal. 
In addition to searching for rare decays, E791 has searched for D° D° mixing ^S], CP violation ["TO] , 
and the pentaquark 20 . We have observed doubly cabibbo suppressed decays [2^, discovered that 
the Df has a longer lifetime than the D° [22], tagged quarks as being c or c using D-ir production 
correlations [23], and measured the — mass splitting [22). Leading effects (e.g. more forward 
ir + (ud) — ► D + {cd) than -K + (ud) — ► D~(cd)) have been observed in D + |2£] and A+ [22j production and 
not in Df production |27| . The total forward and differential cross section for D° production has been 
measured [2Hj. The differential cross section for D* + hadroproduction has been measured [22]. The 
decays A+ -> pK~ir + jHU], D° —* K-K + ir-ir+ [31], D° -> K~K~K + ^ [52], D+ -» K-7r + 7r + [S3], 
and — > 7r~7r + 7r + [SI] have been studied. Evidence for the scalar meson cr(500) in D + — * 7r~7r + 7r + 
decays has been observed [35]. In non-charm physics, E791 has measured A , 3 — , and Q~ hyperon 
production asymmetries jSH) and has used di-jet events to observe the pion valence quark distribution 
|37| and color transparency |38| . 

E791 recorded a total of 2 x 10 10 events with a loose transverse energy requirement and data 
acquisition system writing at 10 MB/s to a great wall of 42 8mm Exabyte tape drives [3JJ|. The 
resulting 50 Terabyte data set was totally unprecedented. It was nevertheless reconstructed at parallel 
computing farms built for this purpose at the University of Mississippi, Kansas State University, 
Fermilab, and CBPF-Rio de Janeiro |40| and yielded an unprecedented 200 000 fully reconstructed 
charmed hadron decays. Track and vertex information came from "hits" in 23 silicon microstrip planes 
and 45 wire chamber planes. This information and the bending provided by two dipole magnets were 
used for momentum analysis of charged particles. Kaon identification was carried out by two multi- 
cell Cerenkov counters [2^ that provided n/K separation in the momentum range 6 — 60 GeV/c. 
We required that the momentum-dependent light yield in the Cerenkov counters be consistent for 
kaon-candidate tracks, except for those in D° — ► (j)ir + ir~ decays with — > K + K~, where the narrow 
mass window for the (ft decay provided sufficient kaon identification (ID). 

Electron ID was based on transverse shower shape plus matching wire chamber tracks to shower 
positions and energies in an electromagnetic calorimeter [22]. The electron ID efficiency varied from 
62% below 9 GeV/c to 45% above 20 GeV/c. The probability to misidentify a pion as an electron was 
~ 0.8%, independent of pion momentum. 

Muon ID employs two planes of scintillation counters [23j- The first plane (5.5 m x 3.0 m) of 
15 counters measured the horizontal position while the second plane (3.0 m x 2.2 m) of 16 counters 
measured the vertical position. There were about 15 interaction lengths of shielding upstream of the 
counters to filter out hadrons. Data from D + — > K \x + v decays [22] were used to choose selection 
criteria for muon candidates. Timing information from the smaller set of muon scintillation counters 
was used to improve the horizontal position resolution. Counter efficiencies, measured using muons 
originating from the primary target, were found to be (99 ±1)% for the smaller counters and (69 ±3)% 
for the larger counters. The probability of misidentifying a pion as a muon decreased with increasing 
momentum, from about 6% at 8 GeV/c to 1.3% above 20 GeV/c. 

Events with evidence of well-separated production (primary) and decay (secondary) vertices were 
selected to separate charm candidates from background. Secondary vertices were required to be 
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separated from the primary vertex by greater than 20 a L for D + decays and greater than 12 a L for D° 
and Dj decays, where a L is the calculated resolution of the measured longitudinal separation. Also, 
the secondary vertex had to be separated from the closest material in the target foils by greater than 
5 a' L , where a' L is the uncertainty in this separation. The vector sum of the momenta from secondary 
vertex tracks was required to pass within 40 /im of the primary vertex in the plane perpendicular 
to the beam. The net momentum of the charm candidate transverse to the line connecting the 
production and decay vertices had to be less than 300 MeV/c for D° candidates, less than 250 MeV/c 
for Df candidates, and less than 200 MeV/c for D + candidates. Finally, decay track candidates were 
required to pass approximately 10 times closer to the secondary vertex than to the primary vertex. 
These selection criteria and kaon identification requirements were the same for both the search mode 
and for its normalization signal (discussed below). 

To determine our selection criteria, we used a blind analysis technique. Before the selection criteria 
were finalized, all events having masses within a window AMg around the mass of the D° were 
"masked" so that the presence or absence of any potential signal candidates would not bias our choice 
of selection criteria. All criteria were then chosen by studying events generated by a Monte Carlo 
(MC) simulation program and background events, outside the signal windows, from real data. 
The criteria were chosen to maximize the ratio Nmc I V Nb , where Nmc an d Nb are the numbers of 
MC and background events, respectively, after all selection criteria were applied. The data within the 
signal windows were unmasked only after this optimization. We used asymmetric windows for the 
decay modes containing electrons to allow for the bremsstrahlung low-energy tail. 

The upper limit for each branching fraction Bx was calculated using the following formula: 



Nx GNorm „ , ^Norm /nottii /,\ 

B X = — T ; x B Novm ; where — = -y^. (1) 

x 



Nx is the 90% confidence level (CL) upper limit on the number of decays for the rare or forbidden decay 
mode X and .BNorm is the normalization mode branching fraction obtained from the Particle Data 
Group • £Norm and ex are the detection efficiencies while and /J^ are the fractions of Monte 

Carlo events that were reconstructed and passed the final selection criteria, for the normalization and 
decay modes, respectively. 

The 90% CL upper limits Nx are calculated using the method of Feldman and Cousins [37] to 
account for background, and then corrected for systematic errors by the method of Cousins and 
Highland |48| . In these methods, the numbers of signal events are determined by simple counting, not 
by a fit. Upper limits are determined using the number of candidate events observed and expected 
number of background events within the signal region. 

Systematic errors include: statistical errors from the fit to the normalization sample iV^orm 

; sta- 
tistical errors on the numbers of Monte Carlo events for both Njf^ and N x lc ; uncertainties in the 
calculation of mis-ID background; and uncertainties in the relative efficiency for each mode, including 
lepton and kaon tagging. These tagging efficiency uncertainties include: 1) the muon counter efficien- 
cies from both Monte Carlo simulation and hardware performance; 2) kaon Cerenkov ID efficiency 
due to differences in kinematics and modeling between data and Monte Carlo simulated events; and 
3) the fraction of signal events (based on simulations) that would remain outside the signal window 
due to bremsstrahlung tails. The large systematic errors for the modes are due to the uncertainty 
in the branching fraction for the Df normalization mode. The sums, taken in quadrature, of these 
systematic errors are listed in Table El and Table 0] 
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The D+ -> h££, D+ Mi, and L>° -> Analysis 

We normalized the sensitivity of our search to the topologically similar Cabibbo-favored decays 
shown in Tableland Figure ^ For the D + decays we used D + — ► K~ir + ir + ; for decays we used 
Df — > i^7r + ; and for D° decays we used D — > i^ _ 7r + events. The efficiencies for the normalization 
modes varied from about 0.5% to 2%. 



Table 1: Normalization modes used for D + — ► -D+ — > and £)° — > £ + i . 



Rare L> Decay 


D Norm. 


Events 


MC Efficiency 


PDG98 46 Branching Ratio 




- 


-* K~TT + TT + 


24010±166 


1.06% 


(9.0 ± 0.6)% 




- 




17730±141 


0.82% 


(9.0 ± 0.6)% 


D+ -> 7r£i 




-> (/>7T + 


952±34 


0.60% 


(3.6 ± 0.9)% 


£>+ -> if+£+r 




-> (/>7T + 


782±30 


0.46% 


(3.6 ± 0.9)% 






-> (/>7T + 


679±27 


0.49% 


(3.6 ± 0.9)% 




L>°- 




25210±179 


1.81% 


(3.85 ± 0.09)% 



The widths of our normalization modes were 10.5 MeV/c 2 for D + , 9.5 MeV/c 2 for D+, and 12 
MeV/c 2 for D°. The signal windows used are: 

1.84 < M(D + ) < 1.90 for D + -> fyi/i 1.78 < M(D + ) < 1.90 GeV/c 2 for D + -► /iee and /i//e 
1.95 < M(£>+) < 1.99 for £>+ -» fyi/i 1.91 < M(D+) < 1.99 GeV/c 2 for D+ -> /iee and fyie 
1.83 < Af(D°) < 1.90 for D° -» ^ 1.76 < M(£>°) < 1.90 GeV/c 2 for £>° -» ee and \ie 





u d) m(K"(lV) GeV/c 2 ° e) m(K~eV) GeV/c 2 ° f) m(K>V) GeV/c 2 

Figure 1: Top row: typical normalization charm signals. The signal region is shaded. Bottom row: invariant 
mass plots of D + candidate decays to K~fi + fi + , K~e + e + , and K~fi + e + , showing reflections mostly from 
misidentified D + — > A"~7r + 7r + decays. These modes are used to set mis-ID rate rather than upper limits The 
solid curves are normalized Monte Carlo fits. The dashed lines show the signal window. 
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Figure 2: Final event samples for the D + (rows 1-3), Df (rows 4-7), and D (row 8) decays. The solid curves 
represent estimated background; the dotted curves represent signal shape for a number of events equal to the 
90% CL upper limit. The dashed vertical lines are AMj boundaries. 

Background that is not removed by cuts include decays in which hadrons (from real, fully-hadronic 
decay vertices) are misidentified as leptons. In the case where kaons are misidentified as leptons, can- 
didates have effective masses which lie outside the signal windows. Most of these originate from 
Cabibbo-favored modes D + — ► K~ir + w + , — ► K~ K + ir + , and D° — ► K~tt + . These Cabibbo- 
favored reflections were explicitly removed prior to cut optimization. There remain two sources of 
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background in our data: hadronic decays with pions misidentified as leptons (A^sId) and "combi- 
natoric" background (Ac m b) arising primarily from false vertices and partially reconstructed charm 
decays. After cuts were applied and the signal windows opened, the number of events within the 
window is N Qhs = ^Sig + ^MisiD + AT Cmb . 



Table 2: E791 90% confidence level (CL) branching fractions (BF) compared to previous experiments. 
The background and candidate events correspond to the signal region only. The Monte Carlo (MC) 
yield is from 250 000 generated events in each of the 24 cases. 







(Est. 


BG) 


Cand. 


Syst. 


90% CL 


MC 


E791 




Previous 




Experi- 


Mode 




A^Cmb 


A^MisID 


Obs. 


Err. 


Num. 


Yield 


BF Limit 


BF Limit 


ment 


D+ -. 


■ ir + fi + fi~ 


1.20 


1.47 


2 


10% 


3.35 


2706 


1.5 x 10" 


-5 


1.8 x 10' 


-5 


E791 P 


D+ -. 


■ 7r + e + e~ 


0.00 


0.90 


1 


12% 


3.53 


816 


5.2 x 10" 


-5 


6.6 x 10 - 


-5 


E791 P 


D+ -. 


- tt + p^e^ 


0.00 


0.78 


1 


11% 


3.64 


1272 


3.4 x 10" 


-5 


1.2 x 1(T 


-4 


E687 [HI 


D+ 


■ 7T~/i + /i + 


0.80 


0.73 


1 


9% 


2.92 


2088 


1.7 x 10" 


-5 


8.7 x 10" 


-5 


E687 [H 


D+ -. 




0.00 


0.45 


2 


12% 


5.60 


701 


9.6 x 10" 


-5 


1.1 x 10- 


-4 


E687 [HI 


D+ -» 




0.00 


0.39 


1 


11% 


4.05 


976 


5.0 x 10" 


-5 


1.1 x 10- 


-4 


E687 [HI 


D + -» 




2.20 


0.20 


3 


8% 


5.07 


1206 


4.4 x 10" 


-0 


9.7 x 10" 


-o 


E687 mi 


D + -» 


K+e + e~ 


0.00 


0.09 


4 


11% 


8.72 


453 


2.0 x 10- 


-4 


2.0 x 10- 


-4 


E687 mi 


-» — V 1 


T r \ _|_ -I- 


0.00 


0.08 


1 


9% 


4.34 


664 


6.8 x 10" 


-5 


1.3 x 10- 


-4 


E687 49 


J*- 


■ K+fi+fi- 


0.67 


1.33 





27% 


1.32 


647 


1.4 x 10- 


-4 


5.9 x 10" 


-4 


E653 E01 




K+e + e~ 


0.00 


0.85 


2 


29% 


5.77 


244 


1.6 x 10- 


-3 












0.40 


0.70 


1 


27% 


3.57 


388 


6.3 x 10- 


-4 












0.40 


0.64 





26% 


1.68 


686 


1.8 x 10- 


-4 


5.9 x 10- 


-4 


E653 [SO] 




K-e+e+ 


0.00 


0.39 





28% 


2.22 


257 


6.3 x 10- 


-4 










K-fi+e + 


0.80 


0.35 


1 


27% 


3.53 


381 


6.8 x 10- 


-4 










7T + fl + ' \l~ 


0.93 


0.72 


1 


27% 


3.02 


1725 


1.4 x 10- 


-4 


4.3 x 10- 


-4 


E653 [SOI 




ir + e + e~ 


0.00 


0.83 





29% 


1.85 


565 


2.7 x 10- 


-4 












0.00 


0.72 


2 


30% 


6.01 


809 


6.1 x 10- 


-4 










TC~ fl + fl + 


0.80 


0.36 





27% 


1.60 


1588 


8.2 x 10- 


-5 


4.3 x 10" 


-4 


E653 [SOI 




ir~e + e + 


0.00 


0.42 


1 


29% 


4.44 


528 


6.9 x 10- 


-4 












0.00 


0.36 


3 


28% 


8.21 


911 


7.3 x 10- 


-4 








D° — > 


/i + fi- 


1.83 


0.63 


2 


6% 


3.51 


5297 


5.2 x 10" 


-6 


4.1 x 10" 


-6 


WA92 [HI 




e + e~ 


1.75 


0.29 





9% 


1.26 


1577 


6.2 x 10- 


-6 


8.2 x 10- 


-6 


E789 mi 


D ^ 




2.63 


0.25 


2 


7% 


3.09 


2983 


8.1 x 10" 


-6 


1.7 x 10" 


-5 


E789 mi 



The background A^MisiD arises mainly from singly-Cabibbo-suppressed (SCS) modes. These misiden- 
tified leptons can come from hadronic shower punchthrough, decays-in-flight, and random overlaps of 
tracks. We do not attempt to establish a limit for D + — ► K~£ + £ + modes, as they have relatively large 
feedthrough signals from copious Cabibbo-favored K~n + ir + decays. Instead, we use the observed 
signals in K~£ + £ + channels to measure three dilepton mis-ID rates under the assumption that the ob- 
served signals (shown in Figs.pi-f) arise entirely from lepton mis-ID. The curve shapes are from Monte 
Carlo. The following mis-ID rates were obtained: r w = (7.3 ± 2.0) x 10~ 4 , = (2.9 ± 1.3) x 10" 4 , 
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and r ee = (3.4 db 1.4) x 10 -4 . Using these rates we estimate the numbers of misidentified candidates, 
^mIsID (f° r D + and Df) and iV^| isID (for D°), in the signal windows as follows: N^f slD = r« • Ngcs 
and N^ islD = r« • N£q S , where iVgQg and iVg^g are the numbers of SCS hadronic decay candidates 
within the signal windows. For modes in which two possible pion combinations can contribute, e.g., 
D + — > h+fj^fjF, we double the rate. 

To estimate the combinatoric background iV Cmb within a signal window AMg, we count events 
having masses within an adjacent background mass window AMb, and scale this number (iV~AM B ) by 
the relative sizes of these windows: iVcmb = (AMj/AMb) • Nam b - To be conservative in calculating 
our 90% confidence level upper limits, we take combinatoric backgrounds to be zero when no events 
are located above the mass windows. 

In Table |21 we present the numbers of combinatoric background, mis-ID background, and observed 
events for all 24 modes. The efficiencies for the search modes varied from about 0.1% to 2%. Data 
are shown in Figure Inland limits are compared with previous results in Figure ^] 

The D° -» V£ + £~ and D° -> hhli Analysis 

There were a few minor differences between this analysis and our previous analysis as discussed 



m 7r+T 



< 150 



above. First, we examined resonant modes, where the mass ranges used were: 

MeV/c 2 , m K ~ 7T + — m-^*o < 55 MeV/c 2 , and \m K + K - — m^\ < 10 MeV/c 2 . We normalized the 
sensitivity of each search to similar hadronic 3-body (resonant) or 4-body (non-resonant) decays. 
One exception is the case of D° — ► p Q l^^ where we normalize to nonresonant D° — > it + it~it + it~ 
because no published branching fraction exists for D° — ► p°-k + it~. Table 01 lists the normalization 
mode used for each signal mode and the fitted numbers of normalization data events (A^Norm)- The 
efficiencies for the normalization modes varied from 0.2% to 1%, and the efficiencies for the search 
modes varied from 0.05% to 0.34%. The signal windows are: 1.83 < M(D°) < 1.90 GeV/c 2 for /jl/j, 
and 1.76 < M(D°) < 1.90 GeV/c 2 for ee and fie modes. 

Table 3: Normalization modes used for D° — > V£ + £~ and D° — > hhil. 



Rare D° Decay D° Norm. Events MC Efficiency PDG2000 53 Branching Ratio 





7T+7T 7T+7T 


2049±53 


0.95% 


(7.3 ±0.5) x 10~ 3 




K TV TV 


5451 ±72 


0.28% 


(1.4 ± 0.4)% 




(pTr + 7r~ 


113±19 


0.21% 


(1.07 ±0.28) x 10~ 3 


ttttU 


7T + 7T~7T + 7T~ 


2049±53 


0.95% 


(7.3 ±0.5) x lO- 3 


KttU 




11550±113 


0.41% 


(7.49 ± 0.31)% 


KKll 


K+K-tt+tt- 


406±41 


0.26% 


(2.5 ±0.23) x 10" 3 



Background sources that are not removed by the selection criteria discussed earlier include de- 
cays in which hadrons (from real, fully-hadronic decay vertices) are misidentified as leptons. These 
misidentified leptons can come from hadronic showers reaching muon counters, decays-in-flight, and 
random overlaps of tracks from otherwise separate decays ( "accidental" sources) . In the case where 
kaons are misidentified as pions or leptons, candidate masses shift below signal windows. However, 
we remove these events to prevent them from influencing our background estimate, which is partially 
obtained from the mass sidebands (see discussion of iVc m b below). To remove these events prior to the 
selection-criteria optimization, we reconstruct all candidates as each of the non-resonant normalization 
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modes and test whether the masses are consistent with m D o . If so, we remove the events, but only if 
the number of kaons in the final state differs from that of the search mode. We do not remove events 
having the same number of kaons, as the loss in acceptance for true signal events would be excessive. 

There remain two sources of background: hadronic decays where pions are misidentified as leptons 
(7V MisID ) and "combinatoric" background (iV"cmb) arising primarily from false vertices and partially 
reconstructed charm decays. The background -/V~MisiD arises from the normalization modes. To estimate 
the rate for misidentifying 7T7t as ££, for all but the D° — > K~ir^£ + £~ modes, we assume all D° — > 
K~tt + £ + £~~ candidates observed (after subtracting combinatoric background estimated from mass 
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Figure 3: Final event samples for the opposite signed dilepton (rows 1-3), resonant (rows 4-6), and same 
signed dilepton modes (rows 7-9) of D° decays. The solid curves display total estimated background; the dotted 
curves display signal shape for a number of events equal to the 90% CL upper limit. The dashed vertical lines 
are the AMs boundaries. 
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sidebands) result from misidentification of D° — > -KT~7r + 7r~7r + decays and count the number of D° — > 
iT — 7T + -£ + i! - decays passing the final selection criteria. We then divide by twice the number of Z?° — ► 
.Rr - 7r + 7r~7r + normalization events with K~ir + £ + l~ mass within AM5 boundaries (twice because there 
are two possible ir + misidentifications). 

Table 0] shows numbers of combinatoric background, misidentification background, and observed 
events for all 27 modes. Data are shown in Figure and compared with previous results in Figure |SJ 

From this procedure, the following misidentification rates were obtained: r w = (3.4 db 2.4) x 10 -4 , 
r fie = (4.2dbl.4) x 10 -4 , and r ee = (9.0±6.2) x 10~ 5 . For modes in which two possible pion combinations 
can contribute, e.g., D° —* if~7r + / u ± /i T , we use twice the above rate; and for D° — > 7r + 7r~7r + 7r~, where 



Table 4: E791 90% confidence level (CL) upper limits on the number of events and branching fraction 
limits (xl0~ 5 ). The Monte Carlo (MC) yield is from 250 000 generated events in each of the 27 cases. 
Previously published limits are listed for comparison |531 1541 150j . 
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Figure 4: Comparison of the 90% CL upper-limit 
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there are 4 possible combinations, we use 4 times this rate in calculating D° —* ir + ir~£ + £~. Using 
these rates, we estimate the numbers of misidentified candidates, N^-J 1D and N^f D , in the signal 
windows as follows: 

A& = r U x ACS and <f s £ ID = r u x < G ™ , (2) 

where and A 7 ^.^ are the numbers of normalization hadronic decay candidates in the signal 

windows. 

To calculate the upper limits for the D° — ► K~ i: + £ + £~ modes, we set A^MisiD to zero as we do not 
have an independent estimate of the misidentification rates. This results in conservative upper limits. 
If we had used the misidentification rates from our previous, 3-body decay study then our limits 
for the three D° — ► K~tt + £ + £~ modes would be lower by about a factor of two. 

To estimate the combinatoric background within a signal window AMg, we count events 

having masses within an adjacent background mass window AMb, and scale this number (N&m b ) by 
the relative sizes of these windows: A/"c m b = (AMg/AMg) x N&m b - To be conservative in calculating 
our 90% confidence level upper limits, we take combinatoric backgrounds to be zero when no events 
are located above the mass windows. 



Conclusion 

We used a blind analysis of data from Fermilab experiment E791 to obtain upper limits on the dilep- 
ton branching fractions for 51 flavor-changing neutral current, lepton- number violating, and lepton- 
family violating decays of D + , Df, and D° mesons. No evidence for any of these 2, 3 and 4-body 
decays was found. Therefore, we presented upper limits on the branching fractions at the 90% con- 
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fidence level. Eighteen limits represented significant improvements over previously published results. 
Twenty-six of the remaining modes had no previously reported limits. Work is currently underway at 
the Fermilab FOCUS experiment and others to further improve the limits presented here or to 
observe signals. 
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